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Abstract 

We present a Wiener filtered map of the Cosmic Microwave Back- 
ground (CMB) fluctuations in a cap with f 5° diameter, centered at the 
North Celestial Pole. The map is based on the 1993-f995 data from 
the Saskatoon (SK) experiment, with an angular resolution around 
1° in the frequency range 27.6-44.1 GHz. The signal-to-noise ratio 
in the map is of order two, and some individual hot and cold spots 
are significant at the 5<r-level. The spatial features are found to be 
consistent from year to year, reenforcing the conclusion that the SK 
results are not dominated by residual atmospheric contamination or 
other non-celestial signals. 



1 Published in ApJ Letters, 474, L77. Available from 
h t t p: //www. mpa-garching.mpg.de/" max/ 'saskmap.html (faster from Europe) 
and from hit p:/ / www.sns.ias.edu/" max/ saskmap.html (faster from the US). 
Note that the figures will print in color if your printer supports it. 



1 INTRODUCTION 



Since the fluctuations in the Cosmic Microwave Background (CMB) depend 
on a large number of cosmological parameters (see Hu, Sugiyama & Silk 
1996 for a recent review), accurate CMB measurements could enable us to 
measure parameters such as the Hubble constant, the density parameter 0, 
etc. to hitherto unprecedented accuracy (Jungman et al. 1996). After the 
successful measurements of large-scale fluctuations by COBE team (Smoot 
et al. 1992; Bennett et al. 1996), attention is now shifting towards measure- 
ments at higher angular resolution. 

When reducing a CMB data set, one usually wants to produce either an 
estimate of the angular power spectrum C\ or a map. Although it is the 
former that is ultimately used to constrain cosmological parameters, there 
are a number of reasons for why map making is useful as well (apart from a 
general desire to map the sky in as many frequency bands as possible): 

• It facilitates comparison with other experiments. 

• It facilitates comparison with foreground templates such as the DIRBE 
maps. 

• It may reveal flaws in the model that are not visible in the power spec- 
trum, such as non-Gaussian CMB features, point sources and spatially 
localized systematic problems. 

The first degree-scale map reconstructed from difference measurements was 
produced by White & Bunn (1995) using data from the MAX experiment, 
which probed a strip of sky a few degrees wide at a resolution of half a 
degree. The purpose of this Letter is to present a map based on data from 
the Saskatoon (SK) experiment (Wollack et al. 1996; Netterfield et al. 1997, 
hereafter "N97"). This map covers a larger patch of sky: a cap of 15° 
diameter centered on the North Celestial Pole (NCP). The angular resolution 
is similar to that of MAX (about half a degree for the SK95 data) , but the 
region in question is more evenly sampled than was the case with the MAX 
map, with no "holes". 

The map-making method we employ is described in Section 2 and the 
results are presented and discussed in Section 3. 

2 METHOD 

The task of generating maps from the Saskatoon data is complicated by the 
fact that the data set does not contain simple sky temperatures, but rather 
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2590 different linear combinations of sky temperatures^. The i th data point, 
which we denote yi, is a linear combination of the temperature across the 
sky, where the weights ascribed to each patch of sky are given by some 
known function /j(r). These weight functions are described in detail in 
N97. For illustration, four sample weight functions are shown in Figure [l]. 
All weight functions emanate approximately radially from the NCP, oscillate 
in the radial direction, cover only a small RA-band, and extend to about 
8° from the pole. The rest of this section describes the inversion process of 
reconstructing a map from these linear combinations 



2.1 Wiener filtering 

Wiener filtering is a general method for estimating a signal from noisy data 
(Wiener 1949), and can be derived as follows.0 Suppose that we have a vector 
of n data points y and wish to estimate a vector of m numbers x (for instance 
the pixel temperatures in a map) from it. Without loss of generality, we can 
assume that both vectors have zero mean, i.e., (x) = and (y) = 0, since 
otherwise, we could redefine them so that they do. Denoting the estimate 
of x by x, the most general linear estimate can clearly be written as 

x = Wy (1) 

for some mxn matrix W. Defining the error vector ciS £ X — X ? cl natural 
measure of the errors is the quantity |e|, which is just m times the r.m.s. 
error per data point. The expectation value of |e| 2 is given by 

(|e| 2 ) = ((Wy-x)'(Wy-x)) = tr [w(yy*)W* - 2(xy*)W* + (xx*)] . (2) 

The Wiener filter W is the matrix that minimizes this error. ^ By dif- 
ferentiating with respect to the components of W, we obtain the simple 
result 

W = (xy i )(yy i )- 1 . (3) 

2 The present analysis is based on the CAP data (as defined in N97) and does not 
include the RING data. 

3 If the data set is Gaussian, an alternative derivation the Wiener filter is to maximize 
the a posteriori probability distribution, i.e., to find the most likely map given the data 
(Zaroubi et al. 1995). 

4 Since the reconstruction of pixel i depends only on the i th row of W, this is equivalent 
to minimizing all the errors e 2 separately. 
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Direct substitution shows that the covariance matrix of the estimates is 

(xx*) = (xy*)(yy*)- 1 (yx i ) (4) 

and that the error covariance matrix is 

(££*) = (xx*)-(xy t )(yy t )- 1 (yx t ). (5) 

Linear filtering techniques have recently been applied to a range of cos- 
mological problems. Rybicki & Press (1992) give a detailed discussion of 
the one-dimensional problem. Lahav et al. (1994), Fisher et al. (1995) and 
Zaroubi et al. (1995) apply Wiener filtering to galaxy surveys. The COBE 
DMR maps have been processed both with Wiener filtering (Bunn et al. 
1994; Bunn et al. 1996) and with other linear filtering techniques (Bond 
1995). 



2.2 Application to the Saskatoon case 

In our case, the observed data point yi is the true sky temperature distri- 
bution x(f) convolved with the i th beam function /«(?), with noise rii added 
afterwards, so 

yi = m + J fi(r)x(r)dtt. (6) 

Our maps will cover a square region of 20° x 20° centered on the NCP, 
pixelized into a 64 x 64 square grid, so m = 4096. To ensure that the maps 
are properly oversampled, we define the pixels to be the sky temperatures 
after Gaussian smoothing on a scale of a = 1°: 

Xi = J <f>(Ti-r)x(r)dn, (7) 
where f j is a unit vector in the direction of the i th pixel and 

4>(cos9)^^e~ e2 ^ 2 . (8) 

This means that the map resolution a is 3.2 times the pixel separation 
20°/64, which is safely above the Shannon oversampling rate of 2.5. In order 
to apply the Wiener filtering procedure, we need to compute the matrices 
(yy*) an d (xy*). The former contains the correlation between the data 
points and themselves, and is given by 

(ViVj) = (niUj) + J J /;(?)/,• (r')c(r • r')dQdn', (9) 



3 



where the correlation function c is given by the angular power spectrum C\ 
through the familiar relation 

00 / 07 _i_ 1 \ 

c(cosfl) = J2 )Pl(cose)C h (10) 

where P; are the Legendre polynomials. As described in N97, the noise 
covariance matrix (nn*) is almost diagonal for the Saskatoon experiment, 
but there are a small number of non-zero correlations (for example, there is a 
2% anticorrelation between the Ka94 5pt East data and the simultaneously 
acquired Ka94 7pt East data due to mainly to atmospheric noise). 

Likewise, the correlation between the data points and the pixels is given 

by 

(x iVj )=J •?')/#)<? -r)dndn'. (11) 

2.3 Practical Issues 

To compute the covariance matrices (xy*) and (yy*), we approximate the 
integrals in equations (^) and ( |TT| ) by sums over a 256 x 265 grid of points. 
Direct computation of (yy*) with this procedure would take over a decade 
on a typical workstation, even if the code were optimized by omitting from 
the double sum all pixels where fi or fj are zero. Fortunately, the relevant 
angular separations are all much less the a radian, which means that the 
effect of sky curvature is negligible. To a good approximation, we can thus 
integrate over a flat two-dimensional plane instead and obtain 

(yy')-(nn')« J J fa) ft (r')c(r - v')d 2 rd 2 r' = J J f t (r)(f^ c){v)d\ 

^ (12) 

where * denotes convolution. Using fast Fourier transforms (FFTs) to com- 
pute the convolutions, the entire covariance matrix can now be computed 
in merely a day. The computation of (xy*) can be accelerated in the same 
way. 

Because the electronic offset is unknown, the means must be removed 
from the observations with each of the 64 synthesized beam patterns, which 
corresponds to multiplying the data vector y by a certain projection matrix 
P. We thus use the corrected covariance matrices (xy*)P* and P(yy*)P* in 
place of (xy*) and (yy 4 ) in equation ([|). We find that this correction makes 
a difference of only a few percent. 
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3 RESULTS AND CONCLUSIONS 



The resulting map is shown in Figure ^ (bottom right) and in Figure ||. The 
fiducial power spectrum used is described in section |3[2| below. As expected, 
it contains virtually no features more than 8° from the center, reflecting the 
fact that the sky outside of this circle was not probed by any of the beam 
functions. Computation of the relevant covariance matrices shows that the 
signal-to-noise ratio within this disc is fairly constant and of order two. In 
other words, the main features visible in this map are expected to be real 
rather than mere noise fluctuations. We also generated a number of mock 
Saskatoon data sets, ran them through the inversion software and compared 
the reconstructions with the original maps, which confirmed this conclusion. 

3.1 Comparison between Years 

The three first panels in Figure |2] show maps generated from the subsets of 
the data that were taken in 1993, 1994 and 1995, respectively. The 1993 data 
is seen to be rather featureless, reflecting the fact that the 1993 data set (42 
data points) contains considerably less information than the other two years 
of data. Similarly, the 1995 map is seen to contain more small-scale struc- 
ture than the 1994 map, reflecting the fact that the angular resolution was 
approximately doubled in 1995. Most potential sources of problems with the 
experiment (underestimation of atmospheric contamination, sidelobe pickup 
from celestial bodies, etc.) would be expected to vary on timescales much 
shorter than a year. In addition, the beam patterns were quite different in 
the three years as described in N97 - for instance, the beam width was sub- 
stantially reduced in 1995 as mentioned above. The visual similarity between 
these independent maps therefore provides reassuring evidence that the bulk 
of the signal being detected is in fact due to temperature-fluctuations on the 
sky rather than unknown systematic problems. 

In addition to a qualitative visual inspection, the 1994 and 1995 maps 
can be used to make more quantitative consistency checks. For example, 
we subtracted one from the other and compared the resulting noise levels 
with the theoretical expectations (as given by the diagonal of W(nn t )W t ). 
The levels are in good agreement, indicating that there is no evidence for 
additional unmodeled/overlooked sources of noise. 

The perhaps most striking single feature in the maps, which stands out 
in all three years of data, is the large cold spot around "two o'clock" , about 
half way from the center. It is interesting to note that the existence of 
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this cold spot can be qualitatively inferred from the plots of the raw data 
in Wollack et al. (1993) and Netterfield et al. (1994), since the three-point 
beam (which has a positive lobe half-way from the center) was found to give 
large negative temperature between 4 and 5 hours in RA. 

3.2 Dependence on method details 

To test if the reconstructed map is sensitive to the pixel size, the analysis was 
repeated with 32 x 32 pixels. As expected, this produced virtually identical 
maps, since the original 64 x 64 pixel map was substantially oversampled. 

The maps in Figure ^ where generated using a a featureless (flat) fiducial 
power spectrum C\ = 6Q 2 /l(l + 1) normalized to Q = 20/iK. To what extent 
do the maps depend on this choice? As described below, the short answer 
to this question is "almost not at all". As a test, we repeated the analysis 
for flat power spectra with Q = 0, 10/uK, 47^K (the best fit to the Saska- 
toon power spectrum points of N97) and 60/iK, as well as for four different 
normalizations of the standard CDM models (Sugiyama 1995). The spatial 
features remained essentially unchanged, and the different normalizations 
simply caused different degrees of smoothing. The appearance of the map 
depended essentially only on one single property of the power spectrum: the 
broad-band power on the angular scales where the SK experiment is sen- 
sitive. This behavior is easy to understand from equation (]3|). Note that 
whereas (yy*) is a sum of two contributions, one from signal and one from 
noise, (xy 4 ) depends only on the signal. Roughly speaking, W is thus of 
the form signal/ (signal+noise). In the extreme case of no signal (Q = 0), 
the Wiener-filtered map thus becomes identically zero, since (xy*) = 0. If 
we increase the assumed signal-to-noise ratio, generic components of W in- 
crease in magnitude, and the Wiener filtering process will attempt to recover 
more details in the map. Since the noise loosely speaking enters on smaller 
scales than the signal, assuming a lower signal-to-noise ratio will basically 
cause the filtering to suppress high frequencies more than low frequencies, 
i.e., smooth the map more. In summary, using a fiducial power spectrum 
with the wrong amount of power in the SK band will produce a map with 
the same spatial features in the same locations, but simply smoothed either 
more or less than what is optimal. 

Which is the best fiducial power level to use? The answer to this question 
depends on our desired signal-to-noise ratio S/N (which we define as the 
ratio of the r.m.s. signal and the r.m.s. noise). The variance in a map pixel 
is given by the corresponding diagonal element of W (yy t )W t , so we can 
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separate the contributions of signal and noise by splitting (yy*) into a signal 
and a noise part and then compute S/N. Since the Wiener filtering balances 
between smoothing too little (getting swamped by noise) and smoothing 
too much (loosing unnecessarily much of the small-scale signal) , it typically 
produces a map where these two problems are comparable in magnitude, 
i.e., where the noise is comparable to the lost part of the signal. Since S/N 
compares the noise to the part of the signal which was not lost, there is 
no a priori guarantee that the S/N obtained will be satisfactory. It is thus 
common to adjust the fiducial power level to obtain a desired S/N. In our 
case, S/N ~ 1.3 for the combined map when the fiducial band power was 
Q = 47{j,K, so we chose Q = 20/uK to get a more smoothed and less noisy 
map, which has S/N ~ 2.0. 

By dividing the map by the r.m.s. noise a, we can read off the significance 
level of individual map features. f\ For instance, the cold spot around "two 
o'clock" is —7(7, the one at "eight o'clock" is —5a and the hot spot at "ten 
o'clock", near the center, is +5cr. 

In conclusion, we have presented the largest map to date of the CMB at 
degree scale angular resolution. The signal-to-noise ratio is of order two, and 
some individual hot and cold spots are significant at the 5a level. It is hoped 
that this map can be used to make comparisons between experiments and 
with various foreground templates, thereby improving our understanding of 
systematics and foregrounds in preparation for the next generation of CMB 
missions. 

The authors wish to thank David Wilkinson for helpful comments on 
the manuscript. Support for this work was provided by NASA through 
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Science Institute, which is operated by AURA, Inc. under NASA contract 
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Forschungsgemeinschaft grant SFB-375, by NSF grant PH 89-21378, NASA 
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5 More generally, equation (|^) can be used to place error bars on linear combinations of 
map pixels (such as correlations with external templates) and to make so-called constrained 
realizations (see e.g. Zaroubi et al. 1995). 
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Figure 1: Sample Saskatoon weight functions. 
Four of the 2590 Saskatoon weight functions are shown in a circle of diameter 
16° with the North Celestial Pole at the center. The weight functions are all 
for the 1995 East data, and correspond to the 3-point beam in the RA =0h 
azimuthal bin (top), the same beam in the RA =6h azimuthal bin (right), 
the 19-point beam in the RA =12h azimuthal bin (bottom) and the 7-point 
beam in the RA =18h azimuthal bin (left). 
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Figure 2: Wiener-filtered maps. 
The CMB temperature is shown in coordinates where the North Celestial 
Pole is at the center of a circle of 16° diameter, with RA being zero at the 
top and increasing clockwise. The first three panels show the maps using 
only the 1993, 1994 and 1995 data sets, respectively. The last panel (bottom 
right) shows the map based on all three years of data. 



Figure 3: The Wiener-filtered map. 
The map based on the entire data set (bottom right in Figure ^) is shown 
in coordinates where the North Celestial Pole is at the center of a 20° x 20° 
square with RA being zero at the top and increasing clockwise. RA = 
is at (x, y) = (0,10°). The contour curves correspond to — 75//K, — 50fiK, 
— 25/xK, 25/iK, 50/iK and 75/xK, respectively. 
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